Molecular characterization and induced changes of histone acetyltransferases in the tick Haemaphysalis longicornis in response to cold stress.

BACKGROUND: Epigenetic modifications of histones play important roles in the response of eukaryotic organisms to environmental stress. However, many histone acetyltransferases (HATs), which are responsible for histone acetylation, and their roles in mediating the tick response to cold stress have yet to be identified. In the present study, HATs were molecularly characterized and their associations with the cold response of the tick Haemaphysalis longicornis explored. METHODS: HATs were characterized by using polymerase chain reaction (PCR) based on published genome sequences, followed by multiple bioinformatic analyses. The differential expression of genes in H. longicornis under different cold treatment conditions was evaluated using reverse transcription quantitative PCR (RT-qPCR). RNA interference was used to explore the association of HATs with the cold response of H. longicornis. RESULTS: Two HAT genes were identified in H. longicornis (Hl), a GCN5-related N-acetyltransferase (henceforth HlGNAT) and a type B histone acetyltransferase (henceforth HlHAT-B), which are respectively 960 base pairs (bp) and 1239 bp in length. Bioinformatics analysis revealed that HlGNAT and HlHAT-B are unstable hydrophilic proteins characterized by the presence of the acetyltransferase 16 domain and Hat1_N domain, respectively. RT-qPCR revealed that the expression of HlGNAT and HlHAT-B decreased after 3 days of cold treatment, but gradually increased with a longer period of cold treatment. The mortality rate following knockdown of HlGNAT or HlHAT-B by RNA interference, which was confirmed by RT-qPCR, significantly increased (P < 0.05) when H. longicornis was treated at the lowest lethal temperature (- 14 °C) for 2 h. CONCLUSIONS: The findings demonstrate that HATs may play a crucial role in the cold response of H. longicornis. Thus further research is warranted to explore the mechanisms underlying the epigenetic regulation of the cold response in ticks.

Unraveling the phylogenetics of genetically closely related species, Haemaphysalis japonica and Haemaphysalis megaspinosa, using entire tick mitogenomes and microbiomes.

Ticks have a profound impact on public health. Haemaphysalis is one of the most widespread genera in Asia, including Japan. The taxonomy and genetic differentiation of Haemaphysalis spp. is challenging. For instance, previous studies struggled to distinguish Haemaphysalis japonica and Haemaphysalis megaspinosa due to the dearth of nucleotide sequence polymorphisms in widely used barcoding genes. The classification of H. japonica japonica and its related sub-species Haemaphysalis japonica douglasi or Haemaphysalis jezoensis is also confused due to their high morphological similarity and a lack of molecular data that support the current classification. We used mitogenomes and microbiomes of H. japonica and H. megaspinosa to gain deeper insights into the phylogenetic relationships and genetic divergence between two species. Phylogenetic analyses of concatenated nucleotide sequences of protein-coding genes and ribosomal DNA genes distinguished H. japonica and H. megaspinosa as monophyletic clades, with further subdivision within the H. japonica clade. The 16S rRNA and NAD5 genes were valuable markers for distinguishing H. japonica and H. megaspinosa. Population genetic structure analyses indicated that genetic variation within populations accounted for a large proportion of the total variation compared to variation between populations. Microbiome analyses revealed differences in alpha and beta diversity between H. japonica and H. megaspinosa: H. japonica had the higher diversity. Coxiella sp., a likely endosymbiont, was found in both Haemaphysalis species. The abundance profiles of likely endosymbionts, pathogens, and commensals differed between H. japonica and H. megaspinosa: H. megaspinosa was more diverse.

Initial study and phylogenetic analysis of hard ticks (Acari: Ixodidae) in Nantong, China along the route of avian migration.

The growing concern about migratory birds potentially spreading ticks due to global warming has become a significant issue. The city of Nantong in this study is situated along the East Asia-Australasian Flyway (EAAF), with numerous wetlands serving as roosting sites for migratory birds. We conducted an investigation of hard ticks and determined the phylogenetic characteristics of tick species in this city. We utilized three different genes for our study: the mitochondrial cytochrome oxidase subunit 1 (COX1) gene, the second internal transcribed spacer (ITS2), and the mitochondrial small subunit rRNA (12 S rRNA) gene. The predominant tick species were Haemaphysalis flava (H. flava) and Haemaphysalis longicornis (H. longicornis). Additionally, specimens of Haemaphysalis campanulata (H. campanulata) and Rhipicephalus sanguineus (R. sanguineus) were collected. The H. flava specimens in this study showed a close genetic relationship with those from inland provinces of China, as well as South Korea and Japan. Furthermore, samples of H. longicornis exhibited a close genetic relationship with those from South Korea, Japan, Australia, and the USA, as well as specific provinces in China. Furthermore, R. sanguineus specimens captured in Nantong showed genetic similarities with specimens from Egypt, Nigeria, and Argentina.

MORPHOLOGICAL AND MOLECULAR CHARACTERIZATION OF AMBLYOMMA SCUTATUM (ACARI: IXODIDAE) ACCIDENTALLY INTRODUCED IN ITALY.

Eight ticks were found in Comacchio (FE), Italy parasitizing a young black iguana (Ctenosaura similis) that had been accidentally transported in a commercial plant container from Costa Rica. Specimens were identified morphologically as Amblyomma scutatum and then confirmed by the barcoding of the mitochondrial cytochrome c oxidase subunit 1 gene. Amblyomma scutatum is a common tick known to infest reptiles in Central America, Mexico, and Venezuela, but not in Europe. In Italy, the possibility for this tick to become endemic is unlikely because of the absence of its principal hosts. Nevertheless, this finding confirms the high risk of introducing exotic species that is linked with global commerce and therefore the need for veterinary control of shipments.

Genetic background of adaptation of Crimean-Congo haemorrhagic fever virus to the different tick hosts.

Crimean-Congo haemorrhagic fever orthonairovirus (CCHFV) is a negative-sense, single-stranded RNA virus with a segmented genome and the causative agent of a severe Crimean-Congo haemorrhagic fever (CCHF) disease. The virus is transmitted mainly by tick species in Hyalomma genus but other ticks such as representatives of genera Dermacentor and Rhipicephalus may also be involved in virus life cycle. To improve our understanding of CCHFV adaptation to its tick species, we compared nucleotide composition and codon usage patterns among the all CCHFV strains i) which sequences and other metadata as locality of collection and date of isolation are available in GenBank and ii) which were isolated from in-field collected tick species. These criteria fulfilled 70 sequences (24 coding for S, 23 for M, and 23 for L segment) of virus isolates originating from different representatives of Hyalomma and Rhipicephalus genera. Phylogenetic analyses confirmed that Hyalomma- and Rhipicephalus-originating CCHFV isolates belong to phylogenetically distinct CCHFV clades. Analyses of nucleotide composition among the Hyalomma- and Rhipicephalus-originating CCHFV isolates also showed significant differences, mainly in nucleotides located at the 3rd codon positions indicating changes in codon usage among these lineages. Analyses of codon adaptation index (CAI), effective number of codons (ENC), and other codon usage statistics revealed significant differences between Hyalomma- and Rhipicephalus-isolated CCHFV strains. Despite both sets of strains displayed a higher adaptation to use codons that are preferred by Hyalomma ticks than Rhipicephalus ticks, there were distinct codon usage preferences observed between the two tick species. These findings suggest that over the course of its long co-evolution with tick vectors, CCHFV has optimized its codon usage to efficiently utilize translational resources of Hyalomma species.

Genotyping of ticks: first molecular report of Hyalomma asiaticum and molecular detection of tick-borne bacteria in ticks and blood from Khyber Pakhtunkhwa, Pakistan.

Multiple ticks (Acari: Ixodoidea) carrying Rickettsiales bacteria have significant importance for both human and animal health. Thus, the purpose of this work was to genetically analyze tick species and their associated Rickettsiales bacteria in animal hosts. In order to achieve these objectives, various animals (including camels, cattle, goats, sheep, dogs, and mice) were inspected in four districts (Mardan, Peshawar, Kohat, and Karak) of Khyber Pakhtunkhwa to collect ticks, while blood samples were collected from all the symptomatic and asymptomatic cattle in all four districts. A total of 234 ticks were obtained from 86 out of 143 (60.14%) host animals, which were morphologically identified as Rhipicephalus turanicus, Rhipicephalus microplus, Haemaphysalis cornupunctata, and Hyalomma asiaticum. Among these, their representative ticks (126/234, 53.85%) were processed for molecular confirmation using cytochrome c oxidase (cox1) gene. Obtained cox1 sequences of four different tick species showed 99.72%-100% maximum identity with their corresponding species reported from Pakistan, China, India, and Kazakhstan and clustered phylogenetically. This study presented the first genetic report of Hy. asiaticum ticks in Pakistan. Moreover, genetically confirmed tick species were molecularly analyzed by PCR for detection of Rickettsiales DNA using partial fragments of 16S rDNA, 190-kDa outer membrane protein A (ompA), and 120-kDa outer membrane protein B (ompB) genes. In addition, blood samples were analyzed to identify Rickettsiales bacteria using the aforementioned genes. Rickettsiales bacteria were found in 24/126 (19.05%) ticks and 4/16 (25.00%) in symptomatic cattle's blood. The obtained ompA and ompB sequences from Hy. asiaticum ticks showed 99.73%-99.87% with Candidatus Rickettsia shennongii and unidentified Rickettsia sp., whereas the obtained 16S rDNA sequences from cattle's blood and ticks (Hae. cornupunctata) showed 99.67% highest identity with Anaplasma phagocytophilum. The 16S rDNA sequence of Rickettsiales DNA from Rh. turanicus ticks showed 100% identity with Ehrlichia canis and unidentified Ehrlichia sp. Obtained sequences of Rickettsiales bacteria were grouped along with their respective species in phylogenetic trees, which were previously reported in Greece, Cuba, Iraq, Turkey, Pakistan, South Korea, and China (mainland and Taiwan). This extensive study explores the wide range of damaging ticks and their corresponding tick-borne bacteria in the area, suggesting a possible danger to both livestock and human communities.

New insights into the molecular phylogeny, biogeographical history, and diversification of Amblyomma ticks (Acari: Ixodidae) based on mitogenomes and nuclear sequences.

BACKGROUND: Amblyomma is the third most diversified genus of Ixodidae that is distributed across the Indomalayan, Afrotropical, Australasian (IAA), Nearctic and Neotropical biogeographic ecoregions, reaching in the Neotropic its highest diversity. There have been hints in previously published phylogenetic trees from mitochondrial genome, nuclear rRNA, from combinations of both and morphology that the Australasian Amblyomma or the Australasian Amblyomma plus the Amblyomma species from the southern cone of South America, might be sister-group to the Amblyomma of the rest of the world. However, a stable phylogenetic framework of Amblyomma for a better understanding of the biogeographic patterns underpinning its diversification is lacking. METHODS: We used genomic techniques to sequence complete and nearly complete mitochondrial genomes -ca. 15 kbp- as well as the nuclear ribosomal cluster -ca. 8 kbp- for 17 Amblyomma ticks in order to study the phylogeny and biogeographic pattern of the genus Amblyomma, with particular emphasis on the Neotropical region. The new genomic information generated here together with genomic information available on 43 ticks (22 other Amblyomma species and 21 other hard ticks-as outgroup-) were used to perform probabilistic methods of phylogenetic and biogeographic inferences and time-tree estimation using biogeographic dates. RESULTS: In the present paper, we present the strongest evidence yet that Australasian Amblyomma may indeed be the sister-group to the Amblyomma of the rest of the world (species that occur mainly in the Neotropical and Afrotropical zoogeographic regions). Our results showed that all Amblyomma subgenera (Cernyomma, Anastosiella, Xiphiastor, Adenopleura, Aponomma and Dermiomma) are not monophyletic, except for Walkeriana and Amblyomma. Likewise, our best biogeographic scenario supports the origin of Amblyomma and its posterior diversification in the southern hemisphere at 47.8 and 36.8 Mya, respectively. This diversification could be associated with the end of the connection of Australasia and Neotropical ecoregions by the Antarctic land bridge. Also, the biogeographic analyses let us see the colonization patterns of some neotropical Amblyomma species to the Nearctic. CONCLUSIONS: We found strong evidence that the main theater of diversification of Amblyomma was the southern hemisphere, potentially driven by the Antarctic Bridge's intermittent connection in the late Eocene. In addition, the subgeneric classification of Amblyomma lacks evolutionary support. Future studies using denser taxonomic sampling may lead to new findings on the phylogenetic relationships and biogeographic history of Amblyomma genus.

Insights from entire mitochondrial genome sequences into the phylogeny of ticks of the genera Haemaphysalis and Archaeocroton with the elevation of the subgenus Alloceraea Schulze, 1919 back to the status of a genus.

We used entire mitochondrial (mt) genome sequences (14.5-15 kbp) to resolve the phylogeny of the four main lineages of the Haematobothrion ticks: Alloceraea, Archaeocroton, Bothriocroton and Haemaphysalis. In our phylogenetic trees, Alloceraea was the sister to Archaeocroton sphenodonti, a tick of an archetypal reptile, the tuatara, from New Zealand, to the exclusion of the rest of the species of Haemaphysalis. The mt genomes of all four of the Alloceraea species that have been sequenced so far had a substantial insert, 132-312 bp, between the tRNA-Glu (E) gene and the nad1 gene in their mt genomes. This insert was not found in any of the other eight subgenera of Haemaphysalis. The mt genomes of 13 species of Haemaphysalis from NCBI GenBank were added to the most recent data set on Haemaphysalis and its close relatives to help resolve the phylogeny of Haemaphysalis, including five new subgenera of Haemaphysalis not previously considered by other authors: Allophysalis (structurally primitive), Aboimisalis (structurally primitive), Herpetobia (structurally intermediate), Ornithophysalis (structurally advanced) and Segalia (structurally advanced). We elevated Alloceraea Schulze, 1919 to the status of genus because Alloceraea Schulze, 1919 is phylogenetically distinct from the other subgenera of Haemaphysalis. Moreover, we propose that the subgenus Allophysalis is the sister to the rest of the Haemaphysalis (14 subgenera) and that the 'structurally primitive' subgenera Hoogstraal and Kim comprise early diverging lineages. Our matrices of the pairwise genetic difference (percent) of mt genomes and partial 16S rRNA sequences indicated that the mt genome sequence of Al. kitaokai (gb# OM368280) may not be Al. kitaokai Hoogstraal, 1969 but rather another species of Alloceraea. In a similar way, the mt genome sequence of H. (Herpetobia) nepalensis Hoogstraal, 1962 (gb# NC_064124) was only 2% genetically different to that of H. (Allophysalis) tibetensis Hoogstraal, 1965 (gb# OM368293): this indicates to us that they are the same species. Alloceraea cretacea may be better placed in a genus other than Alloceraea Schulze, 1919. Reptiles may have been the host to the most recent common ancestor of Archaeocroton and Alloceraea.

Acaricide resistance status of deltamethrin and coumaphos in Hyalomma anatolicum ticks collected from different districts of Haryana.

To study the acaricide resistance status and possible mechanisms of action in conferring resistance to commonly used acaricides (deltamethrin and coumaphos), Hyalomma anatolicum ticks were collected from 6 dairy farms of Hisar and Charkhi Dadri districts of Haryana. By using standard larval packet test, H. anatolicum tick larvae of Charkhi Dadri isolates were found to be susceptible (100% mortality) to both the acaricides. Level-I resistance against coumaphos was recorded from four isolates, whereas, level-II was observed in only one isolate, collected from Hisar. One isolates (Kaimri) from Hisar also showed level-I resistance against deltamethrin. Biochemically, the ticks having higher values of resistance factor (RF) against coumaphos were found to possess increased enzymatic activity of α-esterase, ß-esterase, glutathione-S-transferase (GST) and mono-oxygenase enzymes, whereas, the monoamine oxidase did not show any constant trend. However, the RF showed a statistical significant correlation with GST only. Native PAGE analysis of H. anatolicum ticks revealed the presence of nine types of esterases (EST-1 h to EST-9 h) by using napthyl acetate as substrate. In the inhibitory assay, esterases were found to be inhibited by PMSF, indicating the presence of serine residue at catalytic triad. The partial cds of carboxylesterase and domain II of sodium channel genes were sequenced to determine any proposed mutations in resistant isolates of H. anatolicum ticks, however, no mutations were observed in either gene, indicating that increased expression of detoxification enzymes as a possible mechanism for resistance development, in the current study.

Application of DNA barcodes in the genetic diversity of hard ticks (Acari: Ixodidae) in Kazakhstan.

Forty-five tick species have been recorded in Kazakhstan. However, their genetic diversity and evolutionary relationships, particularly when compared to ticks in neighbouring countries, remain unclear. In the present study, 148 mitochondrial cytochrome c oxidase subunit I (COI) sequence data from our laboratory and NCBI (National Center for Biotechnology Information; https://www.ncbi.nlm.nih.gov/ ) data were used to address this knowledge gap. Phylogenetic analyses showed that i) Hyalomma anatolicum anatolicum (Koch, 1844) ticks from Jambyl Oblast (southeastern Kazakhstan) and Gansu Province (northwestern China) constituted a newly deviated clade; and ii) Dermacentor reticulatus (Fabricius, 1974) ticks from South Kazakhstan Oblast were closer to those in Romania and Turkey. The network diagram of haplotypes showed that i) the H-1 and H-2 haplotypes of Dermacentor marginatus (Sulzer, 1776) ticks from Zhetisu and Almaty were all newly evolved; and ii) the H-3 haplotypes of Haemaphysalis erinacei (Pavesi, 1884) from Almaty Oblast and Xinjiang Uygur Autonomous Region (northwestern China) were evolved from the H-1 haplotype from Italy. In the future, more COI data from different tick species, especially from Kazakhstan and neighbouring countries, should be employed in the field of tick DNA barcoding.

Molecular and next-generation sequencing analysis of tick-borne pathogens of Rhipicephalus ticks (Acari: Ixodidae) in cattle and dogs.

Ticks are small and adaptable arachnid ectoparasites and global carriers of various pathogens that threaten both human and animal health. They are present in many parts of China. A total of 858 ticks were collected from various regions and hosts, then subjected to species identification based on morphological and molecular characteristics, as described in the authors' previous study. Eighty-three individual tick samples were selected for screening pathogens based on metagenomic next-generation sequencing (mNGS) and polymerase chain reaction (PCR) assays. The genomic DNA of tick species was extracted, and amplification of the bacterial 16S rRNA gene was carried out from DNA of individual ticks using V3-V4 hypervariable regions, before subjecting to metagenomic analysis. Each tick underwent specific PCR tests for identifying the bacterial species present, including Anaplasma, Ehrlichia, Coxiella, and Rickettsia, and also protozoans such as Babesia, Theileria, and Hepatozoon. Illumina NovaSeq sequencing results revealed that the dominant phylum and family in Rhipicephalus spp. were Bacteroidota and Muribaculaceae, respectively. Alpha diversity patterns varied depending on tick sex (R. linnaei only), species and location, but not on host. Furthermore, bacterial pathogens, including A. marginale (58 %, 29/50), A. platys (6 %, 3/50), E. minasensis (2 %, 1/50), Ehrlichia sp. (10 %, 5/50), T. sinensis (24 %, 12/50), T. orientalis (54 %, 27/50) and Coxiella-like bacteria (CLB) (80 %, 40/50) were detected in R. microplus, while E. canis (33.33 %, 10/30), H. canis (20 %, 6/30) and CLB (100 %, 30/30) were detected in R. linnaei. Also, Anaplasma sp. (33.33 %, 1/3), A. marginale (33.33 %, 1/3), R. felis (33.33 %, 1/3) and CLB (100 %, 3/3) were detected in R. haemaphysaloides. Dual and triple co-infections involving pathogens or CLB were detected in 84.00 % of R. microplus, 66.66 % of R. haemaphysaloides, and 33.00 % of R. linnaei. The report on microbial communities and pathogens, which found from Rhipicephalus spp. in Hainan Island, is an important step towards a better understanding of tick-borne disease transmission. This is the first report in the area on the presence of Anaplasma sp., A. marginale, R. felis and Coxiella, in R. haemaphysaloides.

Virome diversity shaped by genetic evolution and ecological landscape of Haemaphysalis longicornis.

BACKGROUND: Haemaphysalis longicornis is drawing attentions for its geographic invasion, extending population, and emerging disease threat. However, there are still substantial gaps in our knowledge of viral composition in relation to genetic diversity of H. longicornis and ecological factors, which are important for us to understand interactions between virus and vector, as well as between vector and ecological elements. RESULTS: We conducted the meta-transcriptomic sequencing of 136 pools of H. longicornis and identified 508 RNA viruses of 48 viral species, 22 of which have never been reported. Phylogenetic analysis of mitochondrion sequences divided the ticks into two genetic clades, each of which was geographically clustered and significantly associated with ecological factors, including altitude, precipitation, and normalized difference vegetation index. The two clades showed significant difference in virome diversity and shared about one fifth number of viral species that might have evolved to "generalists." Notably, Bandavirus dabieense, the pathogen of severe fever with thrombocytopenia syndrome was only detected in ticks of clade 1, and half number of clade 2-specific viruses were aquatic-animal-associated. CONCLUSIONS: These findings highlight that the virome diversity is shaped by internal genetic evolution and external ecological landscape of H. longicornis and provide the new foundation for promoting the studies on virus-vector-ecology interaction and eventually for evaluating the risk of H. longicornis for transmitting the viruses to humans and animals. Video Abstract.

Molecular identification of cattle ticks in the Greater Accra Region of Ghana: a high occurrence of Rhipicephalus microplus.

Ticks are competent vectors of a wide range of pathogens. They are of veterinary and public health importance as they affect both animal and human health. Transhumance and the transboundary movements of cattle within the West African Sub-region have facilitated the spread of ticks which threatens the introduction of invasive species. Currently, Rhipicephalus microplus have been identified in the Upper East Region of Ghana which could mean a wider distribution of the species across the country due to livestock trade. This study focused on three sites in the Greater Accra Region, which serves as the gateway to receiving most of the cattle transported from the northern regions of Ghana. Ticks were sampled from August 2022 in the wet season to January 2023 in the dry season. Three tick genera were identified: Amblyomma (19.5%), Hyalomma (1.1%), and Rhipicephalus (79.3%) from the 1,489 feeding ticks collected from cattle. Furthermore, Rhipicephalus microplus, Hyalomma rufipes and Amblyomma variegatum were identified molecularly using primers that target the mitochondrial COI gene. There was a significant association between the tick species and seasons (p < 0.001). Finding R. microplus in this study indicates the extent of the spread of this invasive tick species in Ghana and highlights the need for efficient surveillance systems and control measures within the country.

The diverse functions of Mu-class Glutathione S-transferase HrGSTm1 during the development of Hyalomma rufipes with a focus on the detoxification metabolism of cyhalothrin.

BACKGROUND: Glutathione S-transferases (GSTs) are a superfamily of multifunctional enzymes in living organisms with metabolic and detoxification functions, which can detoxify exogenous and endogenous compounds and thereby reduce the damage caused by toxic substances to the body. Ticks are obligate blood-sucking ectoparasites that can transmit various pathogens, and the characterization of tick-derived GSTs may help improve current understanding of the molecular mechanism of tick resistance to insecticides. In this study, a novel GST gene, named HrGSTm1, was identified from Hyalomma rufipes. METHODS: Sequence analysis was performed by using bioinformatics techniques. A prokaryotic expression system was used to obtain the recombinant expression protein rHrGSTm1. Detection of spatiotemporal expression patterns of target genes and their response to the toxicity of cyhalothrin on female H. rufipes was performed by using a quantitative PCR platform. The optimal enzymological parameters of rHrGSTm1 using glutathione as substrate were calculated. The antioxidant capacity of the recombinant protein was evaluated by DPPH• (1,1-Diphenyl-2-picrylhydrazyl radical 2,2-Diphenyl-1-(2,4,6-trinitrophenyl) hydrazyl). Knockdown of the HrGSTm1 genes through RNA interference was used to analyze their effects on the physiological parameters of ticks. The changes in HrGSTm1 messenger RNA expression patterns under cypermethrin stress were analyzed. RESULTS: The complementary DNA sequence of HrGSTm1 contained a 672-bp open reading frame, which potentially encoded 223 amino acids. The predicted molecular weight was 25.62 kDa, and the isoelectric point 8.22. HrGSTm1 is a Mu-class GST, belonging to the cytoplasmic GSTs with no signal peptide observed. The Vmax and Km of rHrGSTm1 were 3.367 ± 0.81 uM and 2.208 ± 0.76 uM, respectively, and its activities were dependent on different temperatures and pH conditions; the scavenging rate of rHrGSTm1 to DPPH• reached 76.4% at 1.25 mg/ml. Variable expressions of HrGSTm1 were observed under various treatment periods and in different tissues, with the highest appearing in eggs (analysis of variance [ANOVA], F(2, 9) = 279.9, P < 0.0001) and Malpighian tubules (ANOVA, F(3, 12) = 290.5, P < 0.0001). After knockdown of HrGSTm1, compared with the control group, the mortality in the treatment group was increased by 16.7%, the average oviposition rate decreased by 33.9%, the average engorged body weight decreased by 287.38 mg and egg weight decreased by 127.46 mg, although only the engorged body weight was significantly different (t-test, t(44) = 2.886, P = 0.006). After exposure to three sublethal concentrations (LC05, LC10, LC50) of cyhalothrin, the expression level of HrGSTm1 in the midgut, ovary and salivary gland was upregulated, whereas in Malpighian tubules, it showed a trend of upregulation at first and then downregulation, implying different functions during the detoxification in different tissues. CONCLUSIONS: In this study, a novel GST of the Mu-class was successfully isolated from H. rufipes and systematically subjected to bioinformatic analysis and recombination identification. The variation trend of HrGSTm1 expression level in different tissues suggests that the gene has different detoxification functions in different tissues. The potential function of this gene was analyzed to provide basic research for further investigation of its detoxification mechanism.

Comparative analysis of Presence-Absence gene Variations in five hard tick species: impact and functional considerations.

Tick species are vectors of harmful human and animal diseases, and their expansion is raising concerns under the global environmental changes' scenario. Ticks host and transmit bacteria, protozoa and viruses, making the understanding of host-pathogen molecular pathways critical to development of effective disease control strategies. Despite the considerable sizes and repeat contents of tick genomes, individual tick genomics is perhaps the most effective approach to reveal genotypic traits of interest. Presence-Absence gene Variations (PAVs) can contribute to individual differences within species, with dispensable genes carried by subsets of individuals possibly underpinning functional significance at individual or population-levels. We exploited 350 resequencing datasets of Dermacentor silvarum, Haemaphysalis longicornis, Ixodes persulcatus, Rhipicephalus microplus and Rhipicephalus sanguineus hard tick specimens to reveal the extension of PAV and the conservation of dispensable genes among individuals and, comparatively, between species. Overall, we traced 550-3,346 dispensable genes per species and were able to reconstruct 5.3-7 Mb of genomic regions not included in the respective reference genomes, as part of the tick pangenomes. Both dispensable genes and de novo predicted genes indicated that PAVs preferentially impacted mobile genetic elements in these tick species.

Insight into Hyalomma anatolicum biology by comparative genomics analyses.

Hyalomma anatolicum is an obligatory blood-sucking ectoparasite and contributes to the transmission of Crimean-Congo haemorrhagic fever (CCHF) virus, Theileria spp. and Babesia spp. Progress in exploring the adaptive strategy of this ectoparasite and developing tools to fight it has been hindered by the lack of a complete genome. Herein, we assembled the genome using diverse sources of data from multiple sequencing platforms and annotated the 1.96 Gb genome of Hy. anatolicum. Comparative genome analyses and the predicted protein encoding genes reveal unique facets of this genome, including gene family expansion associated with blood feeding and digestion, multi-gene families involved in detoxification, a great number of neuropeptides and corresponding receptors regulating tick growth, development, and reproduction, and glutathione S-transferase genes playing roles in insecticide resistance and detoxification of multiple xenobiotic factors. This high quality reference genome provides fundamental data for obtaining insights into a variety of aspects of tick biology and developing novel strategies to fight notorious tick vectors of human and animal pathogens.

Sequential expression of small heat shock proteins contributing to the cold response of Haemaphysalis longicornis (Acari: Ixodidae).

BACKGROUND: Haemaphysalis longicornis is an important livestock pest and a serious threat to public health. Cold is a common form of stress affecting its survival and distribution. However, H. longicornis exhibits different physiological responses to cold stress. In this study, we systematically explored the regulation and functions of small heat shock proteins (sHsps) in H. longicornis during cold stress. RESULTS: Seven sHsp genes (HlsHsp14.9, HlsHsp19.9, HlsHsp20.3, HlsHsp21.4, HlsHsp23.7, HlsHsp24.0, and HlsHsp26.1) with open reading frame lengths ranging from 408 bp (HlsHsp14.9) to 673 bp (HlsHsp26.1) were cloned from H. longicornis, and featured the typical α-crystallin domain. Phylogenetic analysis revealed high similarity with the sHsps of arachnid species. Quantitative polymerase chain reaction analysis revealed that the regulation of sHsp genes depended on the severity and duration of cold treatment. Moreover, the relative expression of each gene was largely dependent on the treatment period (P < 0.01; 3, 6, and 9 days of treatment at 8, 4, 0, and -4 °C). Among all genes, HlsHsp14.9, HlsHsp19.9, HlsHsp20.3, and HlsHsp24.0 were most sensitive to rapid cold treatment. After RNA interference, the mortality of H. longicornis was significantly increased at -14 °C (P < 0.05), suggesting that the expression of sHsp genes is closely related to cold tolerance in H. longicornis. CONCLUSION: Our results indicate that sHsps play an important role in the cold stress response of H. longicornis, which may enhance our understanding of the cold adaptation mechanisms in ticks. © 2023 Society of Chemical Industry.

Endogenous Viral Elements in Ixodid Tick Genomes.

The documentation of endogenous viral elements (EVEs; virus-derived genetic material integrated into the genome of a nonviral host) has offered insights into how arthropods respond to viral infection via RNA interference pathways. Small non-coding RNAs derived from EVE loci serve to direct RNAi pathways in limiting replication and infection from cognate viruses, thus benefiting the host's fitness and, potentially, vectorial capacity. Here we use informatic approaches to analyze nine available genome sequences of hard ticks (Acari: Ixodidae; Rhipicephalus sanguineus, R. microplus, R. annulatus, Ixodes ricinus, I. persulcatus, I. scapularis, Hyalomma asiaticum, Haemaphysalis longicornis, and Dermacentor silvarum) to identify endogenous viral elements and to illustrate the shared ancestry of all elements identified. Our results highlight a broad diversity of viral taxa as having given rise to 1234 identified EVEs in ticks, with Mononegavirales (specifically Rhabdoviridae) well-represented in this subset of hard ticks. Further investigation revealed extensive adintovirus integrations in several Ixodes species, the prevalence of Bunyavirales EVEs (notably not observed in mosquitoes), and the presence of several elements similar to known emerging human and veterinary pathogens. These results will inform subsequent work on current and past associations with tick species with regard to the viruses from which their "viral fossils" are derived and may serve as a reference for quality control of various tick-omics data that may suffer from misidentification of EVEs as viral genetic material.

Mechanism of the toxic effects of tetracycline on blood meal digestion in Haemaphysalis longicornis.

The extensive utilization of antibiotics in the field of animal husbandry gives rise to various concerns pertaining to the environment and human health. Here, we demonstrate that the administration of tetracycline impedes blood meal digestion in the tick Haemaphysalis longicornis. Tissue sectioning, 16S rRNA high-throughput sequencing, and transcriptome sequencing of the midgut were employed to elucidate the mechanism underlying tetracycline toxicity. The treatment group consisted of engorged female ticks that were subjected to tetracycline microinjections (75 µg per tick), whereas the control group received sterile water injections. On days 2 and 4 following the injections, the tick body weight changes were assessed and the midguts were dissected and processed. Change in tick body weight in tetracycline-treated group was less than in the control group. In tetracycline-treated ticks, midgut epithelial cells were loosely connected and blood meal digestion was impaired compared to the control group. There was no significant change in midgut bacterial diversity after tetracycline treatment. On day 2 following treatment, the relative abundance of Escherichia-Shigella was significantly decreased, whereas the relative abundance of Allorhizobium was significantly increased compared to the control group. On day 4 following treatment, the relative abundance of Escherichia-Shigella, Allorhizobium, Ochrobactrum, and Acidibacter decreased significantly, whereas the relative abundance of Paraburkholderia and Pelomonas increased significantly. Tetracycline treatment also affected midgut gene expression, producing a cumulative effect wherein the differentially expressed genes (DEGs) were mostly down-regulated. KEGG enrichment pathway analysis revealed that on day 2 the up-regulated DEGs were significantly enriched in 21 pathways, including apoptosis and phagosome. Comparatively, the down-regulated DEGs were significantly enriched in 26 pathways, including N-glycan biosynthesis, lysosome, and autophagy. In contrast, on day 4 the up-regulated DEGs were significantly enriched in 10 pathways including aminoacyl-tRNA biosynthesis, ribosome biogenesis, RNA transport, and DNA replication, whereas the down-regulated differential genes were significantly enriched in 11 pathways including lysosome, peroxisome, N-glycan biosynthesis, and fatty acid synthesis. This indicates that tetracycline injection inhibited blood meal digestion by affecting midgut digestive cells, gut flora diversity, and gene expression. These findings could contribute to tick control by inhibiting blood meal digestion.

Transcriptome analysis of Haemaphysalis flava female using Illumina HiSeq 4000 sequencing: de novo assembly, functional annotation and discovery of SSR markers.

BACKGROUND: Ticks are ectoparasites capable of directly damaging their hosts and transmitting vector-borne diseases. The ixodid tick Haemaphysalis flava has a broad distribution that extends from East to South Asia. This tick is a reservoir of severe fever with thrombocytopenia syndrome virus (SFTSV) that causes severe hemorrhagic disease, with cases reported from China, Japan and South Korea. Recently, the distribution of H. flava in South Korea was found to overlap with the occurrence of SFTSV. METHODS: This study was undertaken to discover the molecular resources of H. flava female ticks using the Illumina HiSeq 4000 system, the Trinity de novo sequence assembler and annotation against public databases. The locally curated Protostome database (PANM-DB) was used to screen the putative adaptation-related transcripts classified to gene families, such as angiotensin-converting enzyme, aquaporin, adenylate cyclase, AMP-activated protein kinase, glutamate receptors, heat shock proteins, molecular chaperones, insulin receptor, mitogen-activated protein kinase and solute carrier family proteins. Also, the repeats and simple sequence repeats (SSRs) were screened from the unigenes using RepeatMasker (v4.0.6) and MISA (v1.0) software tools, followed by the designing of SSRs flanking primers using BatchPrimer 3 (v1.0) software. RESULTS: The transcriptome produced a total of 69,822 unigenes, of which 46,175 annotated to the homologous proteins in the PANM-DB. The unigenes were also mapped to the EuKaryotic Orthologous Groups (KOG), Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) specializations. Promiscuous presence of protein kinase, zinc finger (C2H2-type), reverse transcriptase, and RNA recognition motif domains was observed in the unigenes. A total of 3480 SSRs were screened, of which 1907 and 1274 were found as tri- and dinucleotide repeats, respectively. A list of primer sequences flanking the SSR motifs was detailed for validation of polymorphism in H. flava and the related tick species. CONCLUSIONS: The reference transcriptome information on H. flava female ticks will be useful for an enriched understanding of tick biology, its competency to act as a vector and the study of species diversity related to disease transmission.